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ABSTRACT 

We calculate the flux and size ratio of radio emitting features moving in antiparallel 

directions at relativistic speeds, taking into account that the pattern we observe may 

move at a velocity different from the one of the emitting plasma. Both velocities can 

^ ! be determined in sources in which both the approaching and the receding features are 

. observed. This is the case of the galactic superluminal source GRS 1915+105, for which 
O ■ 

the pattern speed (responsible for the apparent superluminal velocity) was found to 

^ I be 0.92c. For this source we find that the velocity of the plasma is 0.73c. The found 

^ ! velocity helps decreasing the estimate of the associated kinetic power. Since intrinsically 

Q_|. identical knots will be observed to have different sizes, high resolution radio observations 

O ■ of GRS 1915+105 will be a test for the proposed difference of the pattern and plasma 

^ ' velocities. 



Subject headings: Relativity - Radio continuum: stars — Stars: GRS 1915+105 



1. INTRODUCTION 



Two oppositely directed blobs of radio emission have been recently observed to move at 
relativistic speeds in GRS 1915+105 (Mirabel & Rodriguez, 1994, hereinafter MR94). 
The apparent velocities are Pa = 1.25+ 0.15 and 0.65+0.08: one of the blobs is therefore 
superluminal. This is the first example that the violent phenomena occuring in active 
galactic nuclei (AGNs) may be reproduced on a much smaller scale also by objects in 
our own galaxy. Due to the proximity of this source, probably located at 12.5 Kpc from 
us, we have a unique oppurtunity to study relativistic effects in much greater depth than 
what is possible for AGNs. Indeed, GRS 1915+105 is the first example of superluminal 
source in which we can observe a pair of blobs going away from a stationary center. This 
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has already allowed MR94 to determine both the real velocity of the moving patterns 
and the angle of motion with respect to the line of sight. Using the derived 'true' speed 
{/3 = 0.92 ±0.08) and viewing angle {9 = 70° ±2°), MR94 were also able to estimate the 
total emitting mass, the kinetic energy of its bulk motion, and the corresponding kinetic 
power, which turned out to be worryingly large: Lf~in ~ lO^"*^ erg s"-*^. Furthermore, the 
knowledge of /? and 6 immediately lead to the prediction of the expected frequency shift 
of any emission line, whose detection would be, of course, very important to confirm 
our ideas about beaming models. 

In this Letter we will argue that the 'true' velocity of the emitting plasma can 
(and very likely is) lower than the velocity of the moving pattern responsible for the 
apparent superluminal speed. This will partly ease the problem in explaining the large 
kinetic power and allows to make a more detailed prediction of the frequency shift of the 
emission line by the outflowing matter. The above results are largely model independent. 
We also calculate the apparent sizes of the approaching and receding knots, and discuss 
how their ratio can give further independent informations on the pattern and plasma 
velocities. 



2.1 FLUX RATIO 

In this section we offer a simple derivation of the expected flux density from two blobs 
moving in antiparallel directions. The pattern of each blob is moving at some velocity 
/3s, while the velocity of the fluid inside each blob is moving with a different velocity Pp. 
All velocities are measured in the observer frame, and the angle 9 is between the line of 
sight and the velocity of the approaching blob. As a special case, we may have Ps = 
(the plasmon case), for which the observed optically thin flux ratio is (for power law 
emission F{i') oc zv"") 



here the subscripts a and r stand for appraching and receding, respectively. 

In the more general case of /3s /3_f, we have to integrate the emissivity (as it is 
in our observer frame) over the observed volume. It is of course better to express the 
emissivity in a frame comoving with the emitting plasma (comoving quantities will be 
denoted by a prime), using 
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where Sp is the Doppler factor of the fluid. A simple way to derive the above transfor- 
mation is to write as the number of photons of energy hu emitted per unit time, 
frequency, solid angle and volume, i.e. as oc dN hu / {dtdiydQclV) and Lorentz transform- 
ing all the quantities (i.e. Rybicki & Lightman, 1979). 

What is the observed volume? When we make a flux measurement (or a map) 
we collect photons arriving at the same time to the detector. For extended objects, 
therefore, these photons leaved the source at different times (Terrel, 1960). For the 
approaching blob, the photons emitted at the back of the blob have to leave at earlier 
times than the photons emitted at the front. This time interval corresponds also to 
different locations of the blob, whose front moves with the pattern velocity /J^c. The 
reverse is true for the receding blob. The ratio of the approaching and receding volumes 
is 

Va ^ 1 + psCOse 

Vr l-PsCOs9 ^ ' 

Therefore the ratio of the observed fluxes is 

Fgjp) ^ l + Ps cose / l + /3FCOs^ y+" 

Special cases are: 

i) (3s = Pf- Then equation (4) reduces to equation (1). 

ii) (3s = 0. This is the case of a stationary jet, whose volume is fixed in the observer 
frame. 

iii) Pf = 0. This is the case of antiparallcl perturbations travelling along the jet, 
powering different portions of plasma at rest (or slowly moving). Even in this case the 
flux ratio is greater than unity, and can reach large values for small angles. 

The same result (i.e. equation 4) can be derived by using the formalism of Lind 
& Blandford (1985), in particular using their equation (2.12b), but we think that our 
derivation is simpler and readily applicable to the case of superluminal sources in which 
both approaching and receding blobs are observed, as is the case of GRS 1915-1-105. 

2.2 SIZE RATIO 

Two opposetely moving blobs, identical in their rest frame, will have different sizes in 
the observer frame. For illustration, consider the case of two cilinders of heigth I and 
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7f(1 ±/3f cos 6) 



(2) 
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section diameter measured in its proper frame, moving at velocities ±/?s along a 
direction which makes an angle d with the line of sight. When we compute its projected 
length seen by the observer, we have to take into account aberration, time delay and 
Lorentz contraction effects. We note that, as in the previous discussion, the projected 
length is that defined by the photons collected at the same time by the observer, which, 
therefore, leaved the source at different times. The total length observed is made by two 
contributions, the projected heigth and the projected section of the cylinder. The first 
one, being parallel to the velocity is modified by Lorentz contraction and time delay 
and becomes 

1°^' = lsme — ^ = ISs sin e (5) 

7s I- Ps cos 6 

For the second one we have only to consider time delay effects, and we have 

d"" = ^^^^^ d (6) 

1 — Ps COS 6 

We note that, when Pg < cos 6, we see the face of the cylinder closer to us, when 
Ps = COS 9, the projected size of the section of the cylinder is zero and we see only its 
lateral surface, while, when Pa > cos 9, we see its farther face. 

The total projected length is 1°''^ + dP^^ and the ratio between that observed for 
the approaching cylinder and that observed for the receding one is 

^ sin6' + X7s|/9s — 0036*1 1 + /3s cos 6* 
sin 9 + x^s (Ps + cos 9) 1 — Ps cos 9 ' 

where x = d/l. The same result could be obtained considering the invariance of the 
perpendicular separation between two parallel light rays (see e.g. the discussion after 
equation 2.6 of Lind & Blandford 1985). Using this method the observed total length 
can be obtained by the following procedure: i) calculate the aberration angle relative 
to the observer and the proper frame; ii) in the proper frame, consider the direction 
of the 'de-aberrated' rays; iii) project the cylinder (in the proper frame) in a plane 
perpendicular to the this direction. This is the projected shape seen by the observer. 
Using this method is immediately apparent that the maximum length of a fast moving 
object cannot be larger than the maximum length seen in the proper frame projecting 
the object in different directions, (i.e. the maximum lenght of a cylinder is Vd? + l^)- 
In a real situation we can imagine that the cylindrical regions observed are asso- 
ciated to shocks and correspond to the regions of emitting post-shock material. Their 
length will therefore be determined by the cooling time of the post-shock material and 
will be in general a function of the frequency of observation. Defining the post-shock 
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plasma velocity, in the shock frame, as and its cooling time, in the plasma proper 
frame, as t^ooli^^)i '^^ compute the cooling length in the shock frame as 

lcool = <ol{^^hFp'F 

In the previous and following expressions the superscripts F, s denote quantities mea- 
sured respectively in the shock and plasma frames. To find how the cooling time de- 
pends on frequency, we must assume a given emission mechanism. For synchrotron 
tcool{j^) which yields t^ooli^^) ~ ^cooli^)^y^^ where f is the observed fre- 

quency. Expressing now Pp and 7!^ in term of the plasma and shock velocities seen by 
the observer and using the observed frequency, we obtain the following expression for 
the quantity x appearing in Eq. (7): 

Special cases are: 

i) (3s — I3f' in this case equation (7) yields R = ^/ {I + Ps)/{^ — Ps) \ cos 9 — 
/3s I /(cos 6* + (3 s), which is always less than unity. 

ii) d <^ ct^ooi '■ ^^^^ limiting case R = {1 + Ps cos 9)/{l — (3s cos 9), indipendently 
on the fluid velocity, reaching very large values for small viewing angles. 

iii) (3f — 0: even in this case we expect R to be, in general, different from unity, 
especially at small viewing angles. 



3. APPLICATION TO GRS 1915+105 

As measured by MR94, the flux ratio of approaching and receding blobs in GRS 
1915+104, calculated at equal separations from the center (hence when the blobs have 
the same age) is ~ 8+1, inconsistent with the value {Fa/Fj. ~ 12 ) expected in the 
case of equal pattern and fluid velocities. Using our equation (4) with Pg = 0.92 and 
9 = 70°, we derive Pf = 0.73 assuming a flux ratio equal to 8 (it becomes 0.664 and 
0.787 with flux ratios equal to 7 and 9, respectively), corresponding to a Lorentz factor 
7F = 1.46. 

Any line emitted by the moving fluid would be shifted in frequency by SF,a = 0.911 
(approaching fluid) and by SF,r = 0.547 (receding fluid), which are signiflcantly different 
by the values found using P = 0.92 and quoted in MR94. 

With the derived value of /3i? we can predict the size ratio R as a function of x{i'), 
which is unfortunately unknown. However, for x{i') =0, 1 and 100 we have R = 0.92, 
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0.84 and 0.68, respectively. Even if the precise value of x{i') is not known, a well defined 
behaviour of R is predicted: R should decrease by observing at different, increasing, 
frequencies, which correspond to shorter cooling times and therefore at larger x{i'). 

Also the estimate of the kinetic power will be significantly affected. In fact, in 
order to estimate the total mass producing the synchrotron emission we observe, we 
have to correct the observed luminosities and frequencies with the correct Doppler 
factor. Following MR94, we may approximate the energy distribution of the emitting 
particles with a monoenergetic function at 7 = 10^73, a tipical value to produce an 
observed frequency of lO^-io^ GHz. This fixes a fiducial value for the magnetic field 
B = 2.7 X 10~^i'iq5/^^. Charge neutrality may be provided by protons or positrons. 
Since electron and positrons would have the same large 7 while protons are very likely 
not relativistic, the hypothesis of a pure pair plasma would only slightly affect (factor 
2) the mass estimate. 

With the values of 5 derived above, we obtain a new estimate for the total mass of 
M = 7.6 X 10^4 g, a factor 3 lower than MR94. The kinetic luminosity is Lkin = 2.6 x 10^° 
erg s~^, a factor ~ 4 lower than MR4. 

DISCUSSION 

Superluminal sources in which both the approcahing and the receding blobs are observed 
offer the oppurtunity to estimate both the velocity of the pattern and the velocity of 
the emitting plasma. For the recently discovered galactic source GRS 1915+105 the 
estimated plasma velocity is Ps = 0.73 ±0.06, significantly lower than the pattern speed 
Ps = 0.92 ±0.08. 

This velocity difference is suggestive of a shock propagating along the jet with 
velocity Ps and postshock velocity Pf- The velocity difference found would indicate 
a compression ratio lower than that expected for strong relativistic shocks in a non 
magnetized medium. Our result is therefore indicative that the magnetic field is likely 
to be important for the dynamics of the flow observed in this source, however more 
precise conclusions are precluded by the uncertainties in the relatistic shock theory 
(Appl & Camenzind 1988; 

If the moving knots are not spheres, as our analysis implies, than their apparent, 
projected shapes are modified by Lorentz contraction and by the different travelling 
times of the photons emitted in different part of the source. This modifications are 
different for the approaching and the receding knots, in such a way that two intrinsically 
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(when at rest) equal knots are observed to have a different total projected size. 

With sufficient angular resolution, the approaching and receding blobs of GRS 
1915+105 should not be observed as circles of equal diameters, but with elongated 
patterns of different maximum extension. The ratio of their sizes can confirm or not 
the plasma velocity derived from the flux ratio. 

As mentioned, any line emitted by the plasma should be observed shifted in fre- 
quency, immediately giving (for a given angle) the corresponding velocity. Formation of 
lines could however be strongly inhibited in a too hot gas. An attractive alternative is 
the annihilation line, expected if part of the jet material is made by electron-positron 
pairs. Note that the source 1E1740-294, close to the galactic center, is thought to 
have jets carrying electron-positron pairs, to explain the variable annihilation line ob- 
served by the SIGMA satellite (Bouchet et al. 1991, Sunyaev et al. 1991, Ramaty et 
al. 1992). GRS 1915+105 may be similar. We may very roughly calculate an upper 
limit to the annihilation line flux assuming that all the particles emitting in the radio, 
once cooled, annihilate before escaping the source. These particles, of initial energy 
7o, emit a fraction (70 — l)/7o of their total energy as radio emission, and a fraction 
1/70 as annihilation radiation. Therefore Fann ~ Fradio/lQ corresponding to a photon 
flux N ~ 10~^^'io-Fjj/(^'io)/7o [photons cm~^ ^~^]- This limit can be increased if only 
a fraction of the existing pairs have been accelerated to high energies to radiate in the 
radio. In this case one can obtain another upper limit by considering that the nucleous 
of the source transforms a fraction ^ of its hard X-ray luminosity into pair rest mass. 
The maximum value of ^ is thought to be around 10 per cent (Svensson, 1987) for non 
thermal plasmas, and much less for thermal ones. In this case, the estimated limit for 
the annihilation photon flux is TV ~ 6.5 x lO^^^Las [photons cm~^ where L38 is 
the luminosity in hard X-rays in units of 10^^ erg s~^. 
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